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Post-Newtonian expansions of the binding energy and gravitational wave flux truncated at the 
same relative post-Newtonian order form the basis of the standard adiabatic approximation to the 
phasing of gravitational waves from inspiralling compact binaries. Viewed in terms of the dynamics 
of the binary, the standard approximation is equivalent to neglecting certain conservative post- 
Newtonian terms in the acceleration. In an earlier work, we had proposed a new complete adiabatic 
approximant constructed from the energy and flux functions. At the leading order it employs the 
2PN energy function rather than the 0PN one in the standard approximation, so that, effectively 
the approximation corresponds to the dynamics where there are no missing post-Newtonian terms 
in the acceleration. In this paper, we compare the overlaps of the standard and complete adiabatic 
templates with the exact waveform in the adiabatic approximation of a test-mass motion in the 
Schwarzschild spacetime, for the VIRGO and the Advanced LIGO noise spectra. It is found that 
the complete adiabatic approximants lead to a remarkable improvement in the effectualness at lower 
PN (< 3PN) orders, while standard approximants of order > 3PN provide a good lower-bound to 
the complete approximants for the construction of effectual templates. Faithfulness of complete 
approximants is better than that of standard approximants except for a few post-Newtonian orders. 
Standard and complete approximants beyond the adiabatic approximation are also studied using 
the Lagrangian templates of Buonanno, Chen and Vallisneri. 

PACS numbers: 04.25Nx, 04.30, 04.80.Nn, 97.60.Jd, 95.55Ym 


I. INTRODUCTION 

Coalescing compact binaries consisting of black holes and/or neutron stars are among 
the most promising sources for ground-based interferometric gravitational wave (GW) 
detectors. One of the main challenges of the data analysis for this kind of sources is to 
create a template bank with which the detector output may be optimally filtered. This, 
in turn, requires an accurate description of the time-evolution of the GW phase. Binary 
coalescences are the end state of a long period of adiabatic dynamics in which the orbital 
frequency of the system changes as a result of gravitational radiation reaction but the 
change in frequency per orbit is negligible compared to the orbital frequency itself. Then 
the inspiral orbit can be thought of as an adiabatic perturbation of a number of circular 
orbits (with a specific conserved energy associated with each of them). Given the binding 
energy E{y) and gravitational wave luminosity fF(v) of the binary, the phasing ip(t) of 
the waves can be calculated in the adiabatic approximation using the following ordinary, 
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coupled differential equations: 

dip 2v 3 dv J~( v ) 

dt m ’ dt mE'(v )’ 

where E'(v ) = dE{v)/dv , and m = mi +m 2 is the total mass of the binary. The binding 
energy -E(u) and gravitational wave luminosity T(v) are calculated, in general, as post- 
Newtonian (PN) expansions in terms of an invariantly defined velocity parameter v . 
The phasing of GWs obtained by numerically solving the above phasing formula is called 
TaylorTl approximant (l|. 


A. Standard and complete approximants 

The standard approach to the GW phasing is based on the PN expansions of the 
binding energy (energy function) and GW luminosity (flux function) truncated at the 
same relative PN order [ 2 ,]. At the lowest order, it uses only the leading terms in the 
energy (Newtonian) and flux (quadrupolar) functions. For higher order phasing, the 
energy and flux functions are retained to the same relative PN orders. We refer to this 
usual physical treatment of the phasing of GWs computed in the adiabatic approximation 
as the standard adiabatic approximation. 

With a view to go beyond the adiabatic approximation, we must think of this in terms 
of the dynamics of the binary under conservative relativistic forces and gravitational 
radiation damping. In the conservative dynamics of the binary, wherein there is no dissi¬ 
pation, the energy is expressed as a post-Newtonian expansion in v 2 , with the dynamics 
involving only even powers of v. When radiation reaction is added to the dynamics, 
then the equation of motion will have terms of both odd and even powers of v. The 
radiation reaction is a correction to the dynamics that first appears at the 2.5PN (i.e. 
v 5 ) order and not at 1PN or 1.5PN order. It is possible to construct the phasing of GWs 
from the direct integration of the equations of motion, without relying on the adiabatic 
approximation. At leading order, the standard non-adiabatic approximant uses the OPN 
and 2.5PN terms in the acceleration (equivalent to using Newtonian conserved energy 
and quadrupolar flux), neglecting the intervening 1PN and 2PN terms. But a complete 
treatment of the acceleration at leading order should include all terms up to 2.5PN, 
without any gaps. We define the phasing of the GWs constructed in this approximation 
as the complete non-adiabatic approximant. 

In the adiabatic approximation, the energy/flux functions can be thought of as car¬ 
rying the information of the conservative/radiative dynamics of the system. So it is 
possible to construct adiabatic variants of standard and complete non-adiabatic approx¬ 
imants entirely based on energy and flux functions. In this sense, standard adiabatic 
approximation is equivalent to neglecting certain conservative terms in the acceleration 
and thus resulting in gaps in the dynamics. In Ref. [|i] , we have proposed a new complete 
adiabatic approximant based on energy and flux functions. At the leading order, it uses 
the 2PN energy function instead of the OPN energy function so that, heuristically, no 
intermediate post-Newtonian terms in the acceleration are missed. Following the nota¬ 
tion used in Ref. Q, we denote the standard adiabatic approximant at nPN order as 
T(£’[ ra ],J r „) and the corresponding complete adiabatic approximant as T(E[ n+ 2 . 5 ],d~n), 
where [p] denotes the integer part of p. 

In Ref. Q, the performance of the standard and complete approximants was evaluated 
in terms of their effectualness and faithfulness 0. We first investigated the problem 
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as a general mathematical question concerning the nature of PN templates assuming 
a flat power spectrum for the detector noise. We then repeated the study using the 
initial LIGO noise spectrum. It was found that, at low (<3PN) PN orders the complete 
approximants bring about a remarkable improvement over the standard approximants for 
the construction of effectual templates. However, the standard approximants are nearly 
as good as the complete approximants at higher orders. In terms of the faithfulness, it 
was again found that the complete approximants generally performed better than the 
standard approximants. 

In the present study, we assess the relative performances of the standard and complete 
adiabatic approximants for the VIRGO and Advanced LIGO noise spectra, as charac¬ 
terized by effectualness and faithfulness. We first perform this for the case of test-mass 
templates by comparing them with an exact waveform calculated in the adiabatic ap¬ 
proximation of a test mass falling in to a Schwarzschild black hole. We then explore 
the extension of the results to the comparable mass case, where the approximants are 
compared with a fiducial exact waveform. We restrict our study to the inspiral phase of 
the coalescing binary, neglecting the plunge and quasi-normal ring down phases. 


B. Noise spectra of the interferometers 


The one-sided noise power spectral density (PSD) of VIRGO is given in terms of a 
dimensionless frequency x = f / fo by Q 


S h (f) = 3.24 x 10 -46 [(6.23a;)- 5 + 2a;- 1 + 1 + x 2 ] 
where fo = 500 Hz; while the same for the Advanced LIGO reads Hi 

1 — x 2 + x a /2 


S h (f) = io - 49 
where fo = 215 Hz. 


x ~4-14 _ 5^-2 _|_ iff ( 


1 + x 2 /2 


( 1 . 2 ) 


(1.3) 


II. TEST MASS WAVEFORMS IN THE ADIABATIC APPROXIMATION 

In the case of a test-particle orbiting a Schwarzschild black hole, the energy function 
E{v) is exactly calculable analytically, while the flux function F(v) is exactly calculable 
numerically |7]. We use these functions to construct the exact waveform in the adiabatic 
approximation. The analytical exact energy function can be Taylor-expanded to get the 
approximants of the energy function. In addition, F(v) has been calculated analytically 
to 5.5PN order jg] by black hole perturbation theory 9]. These approximants of energy 
and flux function are used to construct the approximate templates. In this study, we 
restrict to TaylorTl approximants, since they do not involve any further re-expansion 
in the phasing formula. Hence there is no ambiguity in constructing the phasing of the 
waveforms using approximants with unequal orders of the energy and flux functions as 
required for the complete adiabatic approximant 2 . The exact waveform is also con¬ 
structed by the TaylorTl method. The waveforms (both the exact and approximate) are 
all terminated at tii so = 1/V6, which corresponds to Tl so ~ 86 Hz for the (1M 0 , 50M 0 ) 
binary and F\ so ~ 399 Hz for the (1M 0 , IOMq) binary. The lower frequency cut-off of 
the waveforms is chosen to be F\ ow = 20 Hz. 


2 See Sec. I.B of Ref. jjjj for a detailed discussion. 
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A. Comparison of standard and complete adiabatic approximants 

The effectualness and faithfulness of various PN templates for two archetypical binaries 
with component masses (1Mq, 10M q ) and (1Mq, 50M q ) are tabulated in Tables |T| and 
ITTI All the results are in perfect agreement with the results of our earlier study |.3j 
using the initial LIGO noise spectrum. Complete adiabatic approximants bring about a 
remarkable improvement in the effectualness for all systems at low PN orders (< 3PN). 
On the other hand, the difference in effectualness between the standard and complete 
adiabatic approximants at orders greater than 3PN is very small. Thus, we conclude that 
standard adiabatic approximants of order > 3PN are as good as the complete adiabatic 
approximants for the construction of effectual templates. This study also confirms our 
earlier observation §| that complete adiabatic approximants are generally less ‘biased’ 
in estimating the parameters of the binary. 

Faithfulness of complete adiabatic approximants is generally better at all PN orders 
(even at very high orders) studied, which suggests that the complete approximants are 
closer to the exact solution than the corresponding standard approximants. But there 
are some cases of anomalous behavior. In the next subsection we will try to understand 
the reason for these anomalous cases where the complete approximants perform worse 
than the standard. 


B. Understanding the results 

Table HTil summarizes the PN orders showing the anomalous behavior (i.e. the complete 
approximants being less faithful than the standard approximants) for the different noise 
spectra studied by us. The best-sensitivity bandwidth of each detector is shown in 
brackets 3 . The left-most column in the table shows the flattest noise spectrum and the 
right-most column shows the narrowest one. 

In order to understand the anomalous behavior shown at certain PN orders, we com¬ 
pare the approximants of the tF(v)/E'(v) function with the corresponding exact function. 
Fig. U shows the standard and complete approximants of T(v)/E' (v) along with the ex¬ 
act functions, where a ‘hat’ indicates the corresponding Newton-normalized quantity. 
These results show that, while the complete approximants are far superior to the stan¬ 
dard approximants in modelling the late-inspiral, the early inspiral is better modelled by 
the standard approximants at these PN orders. In the case of the (1M 0 , IOMq) binary, 
the OPN standard approximant is closer to the exact function than the corresponding 
complete approximant in the frequency region 20-50 Hz. But, since none of the detec¬ 
tors is sensitive in this frequency band, this effect shows up in the white-noise case only. 
Similarly the 1.5PN and 2PN standard approximants are closer to the exact function 
in the frequency regions 20-60 Hz and 20-80 Hz, respectively. But the 1.5PN approx¬ 
imant does not show the anomalous behavior in the case of the Advanced LIGO and 
Initial LIGO because the 20-60 Hz region does not fall in the best-sensitivity bandwidth 
of these detectors. The anomalous behavior exhibited by the 3PN approximant can be 
understood in a similar way. It should be noted that the final stages of the inspiral is 
much better modelled by the complete approximants (see the top panel of Fig. [TJ . But, 
since the binary spends more cycles in the early inspiral phase, the overlaps are heavily 
influenced by the efficiency of the modelling of the early inspiral. 


3 It should be noted that there is no rigorous definition for the ‘best-sensitivity’ bandwidth. We define 
it as the bandwidth where the detector’s effective noise amplitude h = fSh(f) is within a factor of 
two of its lowest value. 
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TABLE I: Effectualness of standard(S) adiabatic T(E[ n ],E„) and complete(C) adiabatic 
T{E [n+2 . 5] ,T n ) approximants in the test mass limit. Percentage biases cr m and a v in deter¬ 
mining parameters m and r/ = mirri 2 /m 2 are given in brackets. 


PN Order (n) 


(1M 0 ,1OM 0 


C 


(1Mq, 50M q ) 


c 


Advanced LIGO 




III. COMPARABLE MASS WAVEFORMS IN THE ADIABATIC 

APPROXIMATION 

In the case of comparable mass binaries there is no exact waveform available and the 
best we can do is to compare the performance of the standard adiabatic and complete adi¬ 
abatic templates by studying their overlaps with some plausible fiducial exact waveform. 
The required energy and flux functions have been calculated by supplementing the exact 
functions in the test mass limit by all the known 77 -dependent corrections (from post- 
Newtonian theory ^o|) i n the comparable mass case 4 . In the case of comparable mass 
binaries, the energy function is currentl y k nown u p to 3PN orde r ITU \VA Ibl I hLITt ll?J | 
and the flux function up to 3.5PN order [HE[HMHJ[H|H|3MMMl|^The 
waveforms (‘exact’ and approximate) are constructed by the TaylorTl method and are 
terminated at t;i S o = l/\/ 6 , which corresponds to Pi so ~ 1570 Hz for a (1.4M 0 , 1.4M 0 ) 


4 Sec. IV of Ref. l|| discusses the precise method used in this construction. 
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TABLE II: Faithfulness of standard(S) adiabatic T(Er n ],E n ) and complete(C) adiabatic 
T(E [n+2 . 5] ,T n ) approximants in the test mass limit. 



Advanced LIGO 

VIRGO 


(1M 0 ,1OM 0 ) 

(1M 0 ,5OM 0 ) 

(1M 0 ,1OM 0 ) 

(1M 0 ,5OM 0 ) 

PN Order (n) 

S C 

S C 

S C 

S C 

0PN 

0.1456 0.4915 

0.1608 0.2955 

0.1384 0.3644 

0.1265 0.3881 

1PN 

0.0853 0.1041 

0.1159 0.1609 

0.0682 0.0818 

0.0887 0.1205 

1.5PN 

0.2711 0.3063 

0.2187 0.6735 

0.2524 0.2348 

0.1859 0.5783 

2PN 

0.6998 0.6140 

0.2765 0.8403 

0.7451 0.4617 

0.2514 0.8597 

2.5PN 

0.2143 0.2710 

0.1961 0.3094 

0.2003 0.2496 

0.1612 0.2420 

3PN 

0.8889 0.5791 

0.7252 0.6971 

0.8339 0.5745 

0.7978 0.6210 

3.5PN 

0.7476 0.9985 

0.3852 0.9087 

0.7684 0.9968 

0.3821 0.9259 

4PN 

0.7314 0.8144 

0.4404 0.5761 

0.7501 0.7892 

0.4024 0.5306 

4.5PN 

0.9001 0.9718 

0.5714 0.9078 

0.8753 0.9595 

0.5298 0.9132 

5PN 

0.8273 0.8518 

0.5303 0.6166 

0.8033 0.8232 

0.4968 0.5617 

5.5PN 

0.8376 0.8640 

0.5563 0.6460 

0.8124 0.8340 

0.5147 0.5862 


TABLE III: PN orders showing anomalous behavior in the context of different noise spectra. 
The best-sensitivity bandwidth of each detector is shown in brackets. 


White-noise 

VIRGO 

Advanced LIGO 

Initial LIGO 


(50-400 Hz) 

(60-300 Hz) 

(80-200 Hz) 

(1,1O)M 0 (1,50) 

(1,10)M 0 (1,50)M@ 

(1,1O)M 0 (1, 50)Mq 

(1,10)M 0 (1,50)M 0 

0PN 




1.5PN 

1.5PN 



2PN 

2PN 

2PN 


3PN 3PN 

3PN 3PN 

3PN 3PN 

3PN 


binary and Fi so ~ 220 Hz for a (IOMq, IOMq) binary. The lower frequency cut-off of 
the waveforms is chosen to be Fi ow = 20 Hz. 


A. Comparable mass results in the adiabatic approximation 

Effectualness and faithfulness of various TaylorTl templates in detecting two prototype 
binaries with component masses (10 Mq, IOMq) and (1.4M 0 , 1.4Mq) are shown in Fig. 
13 and are tabulated in Tables ITVl and Ivl It should be noted that the complete adiabatic 
approximants are presently calculated only up to 1PN order. Thus, it is not possible to 
make strong statements of general trends. Heuristically one can conclude from Fig. [3 
that standard adiabatic approximants of order > 1.5PN provide a good lower-bound to 
the complete adiabatic approximants for the construction of both effectual and faithful 
templates. We also note from Table E] that standard adiabatic approximants of order 
2PN/3PN produce the target value 0.965 in effectualness (corresponding to 10% event- 
loss) in the case of the BH-BH/NS-NS binaries. 
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FIG. 1: Top panel shows the approximants of j-(v)/E'(v ) plotted as a function of v. Dashed 
lines indicate standard approximants jF n {y)/E\ n Av) and solid lines indicate the corresponding 
complete approximants F n (v)/E f n+2 5 i(v). Middle and bottom panels show the same approxi¬ 
mants plotted as a function of the GW frequency F = v 3 /irm in the case of the (1 Mg, 10 Mg) 
binary and the (1Mq,50Mq) binary, respectively. 


IV. NON-ADIABATIC APPROXIMANTS 

In we also looked at the relative performance of standard and complete non- 
adiabatic approximants , using the Lagrangian templates described by Buonanno, Chen 
and Vallisneri pflj . comparing them with the exact waveform 5 constructed in the adi¬ 
abatic approximation. The overlaps were computed for the white noise spectrum and 
the initial LIGO noise spectrum. We had found that, as for the adiabatic approximants, 
the complete approximants were far better than the standard approximants for the con¬ 
struction of effectual templates. However, in marked contrast to the adiabatic case, we 
found that the use of complete templates led to a decrease in faithfulness compared to 
the standard templates. 

In this study, we look at the relative performance of the standard and complete ap¬ 
proximants for the case of the Advanced LIGO and VIRGO noise spectra. The results 
are tabulated in Tables l*VTl and I VI J It should be noted that, at present, results are avail¬ 
able at too few PN orders to make statements about general trends in effectualness and 
faithfulness. However, we find that, as in 3|, the complete approximants are generally 
better than the standard approximants for the construction of effectual templates, but 
the faithfulness of the complete approximants is less than that of the standard. 


5 In the comparable mass case, the approximants were compared with a fiducial exact waveform which 
is discussed in Sec. m 


































FIG. 2: Effectualness (four panels in the left) and faithfulness (four panels in the right) of various 
TaylorTl templates in detecting two prototype binaries with component masses (IOMq, 10M@) 
and (1.4 Mq, 1.4Mq). Top panels show overlaps calculated using the Advanced LIGO noise 
spectrum and the bottom panels show overlaps calculated using the VIRGO noise spectrum. 
Different lines in the panels correspond to different orders of the flux function. Each line shows 
how the overlaps are evolving as a function of the accuracy of the energy function. Standard 
adiabatic approximants T(E[ n ]. J- n ) are marked with thick dots. 


V. SUMMARY 

The standard adiabatic approximation to the phasing of GWs from inspiralling com¬ 
pact binaries is based on PN expansions of the binding energy and GW flux truncated at 
the same relative PN order. Viewed in terms of the dynamics of the binary, this standard 
treatment is equivalent to neglecting certain conservative terms in the acceleration. In an 
earlier work pj, we have proposed a new complete adiabatic approximant which, in spirit, 
corresponds to a complete treatment of the acceleration accurate up to the respective PN 
order. In this study we have investigated the performance of the standard and complete 
adiabatic approximants in the cases of the VIRGO and Advanced LIGO noise spectra. 
This has been done by measuring their effectualness (i.e. larger maximum overlap with 
the exact waveform), and faithfulness (i.e. smaller biases in parameter estimation). In 
the test-mass case, the approximants were compared with the exact waveform; while in 
the comparable mass case, the approximants were compared with a fiducial ‘exact’ wave¬ 
form. We have considered only the inspiral phase of the binary, neglecting the plunge 
and quasi-normal ring down phases. The results of this study are in full agreement with 
our earlier studies Q using Initial LIGO noise spectrum and white noise spectrum. We 
summarize the results of our study as follows: 

• In the test-mass case, effectualness of the templates improves significantly in the 
complete adiabatic approximation at lower (< 3PN) PN orders. But standard 
adiabatic approximants of order > 3PN are nearly as good as the complete ap¬ 
proximants. 

• Faithfulness of complete adiabatic approximants is generally better at all PN or¬ 
ders studied. But there are some cases of anomalous behavior. We have shown 
that, at these PN orders (OPN, 1.5PN, 2PN and 3PN) the early inspiral is better 
modelled by the standard approximants than the corresponding complete approx¬ 
imants, which explains the better faithfulness exhibited by the standard approx¬ 
imants at these orders. But complete adiabatic approximants are far superior to 
the standard adiabatic approximants in modelling the final inspiral. 

• Complete adiabatic approximants are generally less ‘biased’ in estimating the pa- 
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TABLE IV: Effectualness of standard(S) adiabatic T(E^, T n ) and complete(C) adiabatic 
T(E\ n+ 2.5lj-7 r n) approximants in the comparable-mass case. Percentage biases a m and a v in 
determining parameters m and r/ = mimz/m 2 are given in brackets. 



(IOMq, 1OM 0 ) 

(1.4M 0 ,1.4M O ) 

PN Order ( n ) 

s c 

S C 

Advanced LIGO 


OPN 

0.7606 (8.5, 0.1) 0.9132 (5.0, 0.3) 

0.8347 (1.4, 0.1) 0.5809 (4.3, 0.1) 

1PN 

0.7110 (57, 0.6) 0.7360 (39, 0.6) 

0.3959 (6.1, 0.0) 0.4194 (5.0, 0.1) 

1.5PN 

0.8741 (2.4, 0.1) 

0.9034 (0.0, 0.2) 

2PN 

0.9803 (0.8, 0.2) 

0.8179 (0.4, 0.0) 

2.5PN 

0.7705 (7.3, 0.1) 

0.7826 (0.4, 0.0) 

3PN 

0.9626 (0.5, 0.0) 

0.9981 (0.4, 0.3) 

3.5PN 

0.9683 (1.3, 1.5) 

0.9977 (0.4, 0.3) 


VIRGO 


OPN 

0.7009 (5.3, 0.0) 0.9280 (4.7, 0.9) 

0.7119 (1.4, 0.1) 0.4405 (5.0, 0.0) 

1PN 

0.5834 (56, 0.3) 0.6148 (30, 0.2) 

0.2808 (3.9, 0.1) 0.2968 (2.9, 0.0) 

1.5PN 

0.8698 (1.3, 0.0) 

0.7724 (0.4, 0.7) 

2PN 

0.9815 (0.8, 0.2) 

0.6420 (0.0, 0.0) 

2.5PN 

0.7299 (4.8, 0.0) 

0.6266 (0.0, 0.1) 

3PN 

0.9624 (0.5, 0.1) 

0.9822 (0.0, 0.3) 

3.5PN 

0.9627 (0.5, 0.1) 

0.9823 (0.0, 0.3) 


TABLE V: Faithfulness of the standard(S) adiabatic T(E[ n ],tF n ) and complete(C) adiabatic 
T(E[ n+ 2 .s], Tn) approximants in the comparable-mass case. 



Advanced LIGO 

VIRGO 


(10Mq, iojv/ 0 ) 

(1.4M 0 , IAMq) 

(IOMq, 1OM 0 ) 

(IAMq, IAMq) 

PN Order ( n ) 

s c 

S C 

s c 

S C 

OPN 

0.3902 0.7030 

0.3731 0.1300 

0.4262 0.5490 

0.3138 0.0794 

1PN 

0.1944 0.2248 

0.1054 0.1128 

0.1574 0.1798 

0.0686 0.0732 

1.5PN 

0.6362 

0.5735 

0.5950 

0.3986 

2PN 

0.8895 

0.3964 

0.8120 

0.3027 

2.5PN 

0.4125 

0.4407 

0.3842 

0.3726 

3PN 

0.9117 

0.9947 

0.9169 

0.9668 

3.5PN 

0.9106 

0.9952 

0.9177 

0.9686 


rameters of the binary. 

• In the case of comparable-mass binaries, standard adiabatic approximants of order 
> 1.5PN provides a good lower-bound to the complete adiabatic approximants for 
the construction of both effectual and faithful templates. 

• In the comparable-mass case, standard adiabatic approximants of order 2PN/3PN 
produce the target value 0.965 in effectualness (corresponding to 10% event-loss) 
in the case of the BH-BH/NS-NS binaries. 
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TABLE VI: Effectualness and faithfulness of the standard(S) and complete(C) non-adiabatic 
approximants in the test mass case. 



Effectualness 

Faithfulness 


(1 Mq, 10 Mq) 

(1M 0 ,5OM O ) 

(1 M @ , 10Mq) 

(1M 0 ,5OM 0 ) 

PN Order (n) 

s c 

S C 

S C 

S C 

Advanced LIGO 

0PN 

0.4259 0.8682 

0.6384 0.9360 

0.2636 0.0754 

0.4857 0.0999 

1PN 

0.5256 0.8280 

0.6080 0.9211 

0.3419 0.1368 

0.4493 0.2038 

VIRGO 

0PN 

0.3720 0.7631 

0.5985 0.9618 

0.1991 0.0570 

0.4947 0.1057 

1PN 

0.3599 0.7386 

0.5777 0.9525 

0.2499 0.0911 

0.5259 0.1954 


TABLE VII: Effectualness and faithfulness of the standard(S) and complete(C) non-adiabatic 
approximants in the comparable mass case. 



Effectualness 

Faithfulness 


(1OM 0 ,1OM 0 ) 

(1.4M 0 ,1.4M 0 ) 

(1OM 0 ,1OM 0 ) 

(1.4M 0 ,1.4M 0 ) 

PN Order (n) 

S C 

S C 

S C 

S C 


Advanced LIGO 


0PN 

0.9147 

0.4338 

0.7417 

0.8322 

0.0637 

0.0546 

0.4662 

0.2192 

1PN 

0.3937 

0.5132 

0.7443 

0.8158 

0.0794 

0.0703 

0.6594 

0.4788 

VIRGO 

0PN 

0.8142 

0.3113 

0.6895 

0.7341 

0.0414 

0.0334 

0.3439 

0.1661 

1PN 

0.2880 

0.3944 

0.6807 

0.7420 

0.0803 

0.0463 

0.5709 

0.3704 
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